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On a Class of Spatial Discretizations of Equations of the Nonlinear Schrodinger Type 
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We demonstrate the systematic derivation of a class of discretizations of nonlinear Schrodinger 
(NLS) equations for general polynomial nonlinearity whose stationary solutions can be found from 
a reduced two-point algebraic condition. We then focus on the cubic problem and illustrate how 
our class of models compares with the well-known discretizations such as the standard discrete NLS 
equation, or the integrable variant thereof. We also discuss the conservation laws of the derived 
generalizations of the cubic case, such as the lattice momentum or mass and the connection with 
their corresponding continuum siblings. (This manuscript was submitted for publication on October 
14, 2005.) 
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I. INTRODUCTION 

In the past few years, the role of spatial discreteness in 
lattice systems described by differential-difference equa- 
tions has been increasingly recognized Q. In these set- 
tings, the spatial variables are discrete, while the evo- 
lution variable is continuum. Relevant applications are 
continuously arising in rather diverse physical contexts 
such as the spatial dynamics of optical beams in coupled 
waveguide arrays in nonlinear optics 0, the temporal 
evolution of Bose-Einstein condensates (BECs) in opti- 
cal lattices in soft-condensed matter physics 3J , the DNA 
double strand in biophysics 4], and so on. 

In the examples that stem from physical applications, 
the form of the discrete model is dictated by the under- 
lying physics, and typically that form is the discrete non- 
linear Schrodinger (DNLS) equation [s^. However, there 
are also exceptions to this rule, where the nature of the 
nonlinearity Q or of the dispersion or both Q imposes 
variations of this ubiquitous model. Another motivation 
to study such modified DNLS models is a more mathe- 
matical one, namely the aim of identifying models with 
good mathematical properties (e.g. exact solutions, ad- 
ditional symmetries or, possibly, complete integrability) . 
Such a program was initiated by the derivation of an 
integrable analog of the DNLS equation, namely the so- 
called Ablowitz-Ladik model |^. This was later used for 
computational studies of the NLS equation T^, as well 
as implemented as a good starting point for developing 
perturbation theoretic approaches to the DNLS limit, in 
order to examine the existence and stability of its solu- 
tions 11] . 

We should note here that a similar, mathematically- 
minded program of discretizations has been evolving 
for Klein-Gordon (KG) type lattices. For instance, in- 
tegrable discretizations of equations such as the sine- 
Gordon model have been obtained [13. H^ : there have 



also been attempts to systematically discretize preserv- 
ing symmetries of the underlying continuum model, such 
as, especially, a discrete analog of translational invari- 
ance. It has often been noted that "standard" (e.g. 
centered-difference schemes for Laplacian type operators) 
discretizations strongly violate translational invariance, 
leading stable and unstable steady states (typically cen- 
tered on-site and inter-site between two lattice nodes). 
This not only imposes undesirable mathematical proper- 
ties |T^ . but also modifies the underlying phenomenol- 
ogy in comparison to the continuum model |15i | . In that 
view, discretizations that preserve an effective transla- 
tional invariance by allowing the center of a stationary 
state to be a free parameter (rather than to be fixed on- 
or half-way between two lattice sites) have been sought. 
In this way, such discretizations also avoid "energy barri- 
ers" (so called Peierls-Nabarro barriers) between on-site 
and inter-site states. Such discrete models have been 
constructed in the KG case, based on a discretization 
of the energy, using the Bogomol'nyi approach 0, as 
well as a discretization of the equation of motion ensur- 
ing the persistence of a discrete momentum conservation 
law 0, . These classes of models were subsequently 
tested for the potential bearing of travelling wave solu- 
tions, using the technology based on the calculation of 
the Stokes constant 193; this led to the conjecture that 
such models may possess isolated, exact, travelling lat- 
tice solutions. It is, finally, worthwhile to note that while 
motivated by their mathematical properties, such models 
may also bear physical relevance as is indicated e.g. in 
the very recent preprint of |23| for the sine-Gordon case. 

While this technology has been well developed for the 
single (i.e., scalar) field case of KG lattices, such con- 
siderations do not seem to have been applied to NLS 
type lattices, to the best of our knowledge. The present 
manuscript aims to partially fill this gap, by present- 
ing a systematic methodology for deriving discretizations 
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of polynomial nonlinearity partial differential equations 
(PDEs) of the NLS type. The main novel feature of these 
discrete models is that, contrary to what is the case for 
the standard DNLS equation, they preserve a discrete 
analog of the momentum conservation law. In fact, we 
show that in the cubic case, they are natural homotopic 
generalizations of the integrable NLS discretization of Q ■ 
Our presentation will be structured as follows. In sec- 
tion II, we present the general setup of the continuum and 
discrete models of the present work. In section III, we 
illustrate the auxiliary problem that aids us to construct 
the desired discretizations in section IV. In section V, 
we study some of the conservation laws of the obtained 
models, while in section VI, we discuss their solitonic 
properties. Finally, in section VII, we briefly summarize 
our findings and present our conclusions. 



II. SETUP 

We present our methodology for the generalized NLS 
equation of the form 



(1) 



where ip{x, t) is a complex function of two real variables; 
G(^) is a real function of its argument and G'{£,) = 
dG/d^. 

We introduce the lattice x„ — nh, where h is the lat- 
tice spacing and n = 0, ±1, ±2, ... We also introduce the 
following shorthand notations 



V'„-i = and Vn+i = V'-H, 



and will focus only on discretizations that involve such 
nearest neighbor sites. 

Our more specific aim will be to construct the discrete 
analogues of Eq. ^ of the form: 

i'(/'„ -Hr('0-,V'-,V'n,V'n,'0+,V'+) = 0, (3) 

such that the ansatz 

Mt) = /™e^"*, (4) 

reduces Eq. (|3J to the three-point discrete problem of 
the form 



^/«+i?(/-,/n,/+) =0, 



(5) 



whose solution can be found from a reduced two-point 
discrete problem m(/_,/„) = 0. Such a selection will 
entail a mono-parametric freedom for the resulting alge- 
braic problem leading to stationary state solutions. This 
will, in turn, be responsible for the effective translational 
invariance in what follows. 



III. AUXILIARY PROBLEM 

Firstly, we formulate an auxiliary problem. Seeking 
stationary solutions of Eq. JQ) in the form 



(6) 



we reduce it to an ordinary differential equation (ODE) 
for the real function f{x), 

D{x) = }" -2uf + 2fG'{f) = 0, (7) 

having the first integral 

u{x) EE (/')' - 2ujf + 2G{f) = 0. (8) 

We then identify discretizations of Eq. ([7|l of the form 

^(/-,/n,/+)=0, (9) 

such that solutions to the three-point discrete Eq. 
can be found from a reduced two-point problem 



■U{f-,fn) = Un - f-f 

-2c./_/„ + 2G(/2,/2) = 0, 



(10) 



which is a discrete version of Eq. 0, assuming that 
G{fl, fn) reduces to G{p) in the continuum limit {h 
0). 

Taking into account that Eq. {Tj) is the static Klein- 
Gordon equation with the potential 



n/) = ^/'-G(/2), 



(11) 



(^2) a wide class of discretizations solving the auxiliary prob- 
lem has been offered in the very recent work of |18l |. 

For example, discretizing the left-hand side of the iden- 
tity {l/2)du/df — D{x), we obtain the discrete version 
ofEq. CZI), 



Jn, J + ) = 7 7 = U. 



(12) 



Formally, Di{f-, fn, f+) = is a three-point problem 
but, clearly, its solutions can be found from the two-point 
problem u{f^, /„) = and thus, the auxiliary problem is 
solved. We note, in passing, that this type of argument 
was first proposed in 



IV. MAIN PROBLEM 

Coming back to our main problem of finding special 
discretizations for Eq. |^ , we should remark that among 
the solutions to the auxiliary problem we should select 
the ones which can be rewritten in terms of ipn and ip^ 
the desired form of Eq. This can be done when Di 
given by Eq. H12|) is written in a non-singular form (i.e., if 
the denominator cancels with an appropriate factoring of 
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the numerator). This always occurs if G(^) is polynomial 
and if u{f_, /„) possesses the symmetry 

"(/-,/«)=«(/«,/-)• (13) 
We thus focus on G(^) in the form of Taylor expansion, 



(14) 



k=l 



with real coefficients a^; retaining first four terms of the 
expansion, we write 



(/l + /,t) + (l-")/^/^ 



+02 



+ 04 



+ (l-7-<5)/l/„* 



(15) 



where a, (3, 7, and 5 are free parameters. The symmetry 
condition of Eq. (|13|l is satisfied for u{f^,fn) given by 
(Unj since in Eq. (O we have G{flJ^) = G{fl, fi). 

Equation ifT^ with «(/-,/«) given by Eq. ((Tn|l and 
G{fl, /^) given by Eq. (|15() assumes the following form 



-/„ + 2p (/- 



2/n + /+) 
+ (/_+/+)i?(/!,/2,/2)=0. 



(16) 



where 



ai 



0.2 



if- + fl) + (1 - «)/^ 



y/3(/l+/^/U/+) 



+ y(l-mU/-+/« + /^) 
7 
2 
(5 



+04 (/!+/:) [j {n+n) + {i-i-5)fn 



One can conclude that the discretization of NLS equa- 
tion in the form: 

1 



(18) 



with R given by the expression of H17|l satisfies the gen- 
eralized equation |^ with G{£) given by Eq. (|14|l . How- 
ever, as per the construction above, additionally, the cor- 
responding discrete equation for the stationary solutions 
of the form of Eq. is the three-point problem of Eq. 
(|16|l , whose solution can be found through the two-point 
reduction of Eq. itTUI) where G(/^,/^) is given by Eq. 



It is interesting to note that the "standard" DNLS 
equation 



2i/;„ + i/;+) + |7/>„|Vn = 0, (19) 



does not belong to the above class and more generally 
does not share the reduction property used above. In- 
stead, and focusing only on the cubic Kerr nonlinearity 
(where G' is linear in its argument), we obtain from Eq. 
(|T5|l and Eq. in|) 



1^4 



(1 - a)\'ip„ 



(20) 



Notice that the integrable discretization of is obtained 
from this approach as the special case of a = 0. For 
a 7^ 0, this model can be regarded as a Salerno- type 
model j2ll |. i.e., a homotopic continuation including the 
integrable limit and reducing to NLS equation in the con- 
tinuum limit. 

Expression (|15|l contains only the terms with even pow- 
ers of /_ and /„. However, it is possible to construct the 
terms of desired symmetry involving odd powers |lS| . For 
example, for Kerr nonlinearity one can take 



G(/-,/„)--/-/„(/! + /^), 
and obtain from Eq. (^HJ and Eq. ifT^ 



- ^/n + ^ (/- - 2/„ + /+) 



(21) 



4 + T (/- + /-/+ + /+) 



0, 



(22) 



for which the following discretization can be obtained 

+ i(|^_|2 + |V.-V+l + l^+l')^n = 0. (23) 

The model of Eq. (|23|) has on-site cubic nonlinear- 
ity modified through inter-site coupling, which makes it 
qualitatively different from integrable system of . Soli- 
ton solutions to this model can be constructed from the 
quartic Eq. ((Tn|l with G given by Eq. (|^ . Note that 
DNLS equation with the quintic term possessing a struc- 
ture similar to the last term of Eq. H23|) has been consid- 
ered in [see Eq. (74) of that paper]. 



V. MOMENTUM AND MASS CONSERVATION 
LAWS 

We now try to connect the above presented construc- 
tion to the relevant conservation laws of the resulting 
infinite-dimensional dynamical system. More specifically. 
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we examine the momentum conservation law that, as ar- 
gued above, is intimately related to the translational in- 
variance and the existence of mono-parametric stationary 
solutions. For the DNLS model of Eq. we consider 
the momentum defined as 



n— — oo 

oo 



(24) 



n— — oo 



We now demand that the momentum be conserved, 
i.e., that 



dP 

dt 



= 0. 



(25) 



Upon substitution of Eq. (|24() and use of the equation of 
motion Eq. i-^n -f r„ — 0, for Vjn, V'^ and -01, we 
obtain 

oo oo 
n— — oo n— — oo 

In the last expression, if the first sum is zero, then the 
second sum is also zero, being its complex conjugate. 
Thus, a sufficient condition for the conservation of the 
momentum is 



(27) 



A direct consequence of the above is that the regular 
DNLS equation, will not succeed in leading to conserva- 
tion of momentum, while it can easily be checked that 
the opposite is true for the model of Eq. 118|l with R 
given by Eq. ((T7|) . 

The Eq. (|18|) with R given by Eq. H17|) conserves mo- 
mentum Eq. H24|) but it does not conserve the "standard" 
(P) norm 



N 



(28) 



In fact, it can be shown that unless R is "local" (i.e., de- 
pendent on l-iAnP and not its neighbors, as is e.g. the case 
for the integrable discretization of ) , there is no defini- 
tion of = X]n-^(I^"P) which can be preserved under 
(|18|l . Instead Eq. preserves the "mass" (norm) of 

the form: 



(29) 



which in the continuum limit retrieves the standard con- 
servation law of the norm. Notice, however, that the 
discretization of Eq. H23|) does conserve the standard 
norm. 



For reasons of comparison that will become more trans- 
parent below, let us also introduce an additional dis- 
cretization that does not belong to the family of Eq. H18|) 
presented above: 



1 



1 



+-^{r-+r+X = o. (30) 



The particular feature of this dynamical system is that it 
conserves the norm but does not conserve the momen- 
tum defined above. However, considering the standing 
wave ansatz of Eq. Q reduces Eq. to 



1 

2^ 



(/- - 2/„ + f+) + \ if- 



+ /+) fn = 0(31) 



which is particular case of Eq. H16|l at 02 = 1/2, ai = 
03 = — and a = 0. Thus, solution of Eq. (|31|l can 
be found from the two-point difference equation 



U{f-Jn) = {fn 



-2LuUU + fif^^0, (32) 



which is particular case of Eq. pOfl . As a result, 
the stationary states can be found by the solution of 
a reduced two-point problem; namely, for any /_ (or 
fn) i n the range where qm = V^io and Qs = 

h^^^/T~~{T+TJhFy^, solving the binomial Eq. Ipl^ . 
one can find /„ (or /_), thus reconstructing the soliton 
solution for Eq. H30|l in the form of Eq. Q). Quantities 
Qm and Qs are the amplitudes of solitons centered be- 
tween two lattice sites and on a lattice site, respectively. 
Interestingly, qm is not a function of h. Even more inter- 
estingly perhaps, such stationary solutions of Eq. 
are identical (see also below) to those of the integrable 
equation [Eq. (^0)1 at a = 0] since the two models share 
the same reduced two-point problem [Eq. 



VI. COMPARISON OF SOLITON SOLUTIONS 
FOR DIFFERENT DISCRETIZATIONS 

We now compare some properties of the classical DNLS 
model of Eq. ifT^ll (model /), the "secular" model of 
Eq. H30() conserving the classical norm (model //), and 
the Kerr-representativc of the class of models developed 
herein, given by Eq. H2(JI) (model III). 

All three models share the same continuum limit, the 
integrable NLS equation with Kerr nonlinearity. 



1 



(33) 



and thus, in the regime of weak discreteness (small lattice 
spacing /i), their soliton solutions of the form of Eq. Q 
can be expressed approximately as 



cosh[qh{n — xq)] 



exp[-z(9V2)i], (34) 



5 



where q and lo = are the sohton amphtude and 

frequency, respectively. 

The approximate solution of Eq. H34|l contains the free 
parameter xq defining the soliton position. However, as 
indicated above, in contrast to the NLS equation of Eq. 

where a;o can be chosen arbitrarily due to transla- 
tional invariance, the DNLS models usually have station- 
ary soliton solutions only for a discrete set of values of xq 
(e.g. on-site and inter-site, as mentioned above). This is 
true, for example, for the classical DNLS of model / and 
for the Salerno model among others. The models 
// and by construction, are among the members of 
a wider class of DNLS equations proposed in this paper, 
where stationary soliton solutions exist for any xq, or, in 
other words, they can be placed anywhere with respect to 
the lattice; otherwise put, the Peierls-Nabarro potential 
is absent for stationary solutions of these models. 

An explicit formula, as is well-known does not ex- 
ist for the solutions of model /. Such solutions can be 
obtained numerically with the desired degree of accu- 
racy, for two particular cases of = and Xq = 1/2 
(due to the integer shift-invariance of the lattice, we now 
restrict ourselves to < < 1). To obtain the soli- 
ton of frequency lu centered on a lattice site, i.e., the 
one with xq — 0, we set fo = \p2xo which is the soli- 
ton amplitude estimated from Eq. (|34|) . Having /o and 
the symmetry property /„ = /_„ for n > 0, we find 
successively /i, then /2 and so on from the equation 
-w/„ + (2/i2)-i (/_ _ 2/„ -f /+) + /3 = 0, which is ob- 
tained by substituting Eq. Q into Eq. (fT^ . Since /q is 
not an exact value of the soliton amplitude, the boundary 
conditions /_oo = /oo = will not be satisfied. We find 
numerically a correction to /o to satisfy the boundary 
conditions thus completing the construction of the soliton 
solution. For xq different from or 1/2 it is impossible 
to satisfy both boundary conditions simultaneously. The 
soliton with xq = 1/2 can be constructed similarly us- 
ing the symmetry property /_i = /o and /„ = /-„-i for 
n > 0. Estimation of the soliton amplitude from Eq. H34|) 
in this case is /o = \/2uj/ cosh[(/i/2)\/2u;]. Alternatively, 
one can straightforwardly use fixed point algorithms to 
obtain such solutions as is summarized in Q . 

For the models // and III the exact solutions of the 
form of Eq. Q can be found using the method developed 
in section IV. More specifically, however, for model //, 
as is expected from the discussion above and the coinci- 
dence of the reduced two-point problem with that of the 
integrable discrete model, an exact stationary solution 
can be obtained explicitly in the form 



1 



sinh /X 



h cosh[/x(n — xq)] 



exp" 



(35) 



where xq is the parameter defining the soliton posi- 
tion and it can obtain any value from [0, 1). The soli- 
ton frequency uj = h~^{l — cosh/i) and amplitude q — 
h^^ sinh^ ^ are expressed in terms of the free parameter 
/i > 0. The surprising feature of this secular discretiza- 
tion is that, despite the absence of an explicit momentum 




FIG. 1: The dynamical evolution of IV'nP for exact solitary 
wave solutions in the classical DNLS model I centered at (a) 
Xq = and (b) xq = 1/2. Soliton in (a) is stable while in (b) 
it is unstable and, due to the presence of perturbations in the 
form of round-off errors, it spontaneously starts to alternate 
between two nearest inter-site configurations passing through 
the stable on-site configuration. Results for lattice spacing 
h = 0.4 and soliton frequency u — 1. 



conservation law, the stationary solutions appear to en- 
joy an effective invariance with respect to their center 
location. 

Model /// has the solutions of the form of Eq. Q 
with /„ derivable from the two-point problem 



1 



(/- - fnf 



2c./_/„ + ^(/l 



/,t) 



+ (1 - a)/2 /2 = 0. 



(36) 



Equation (|36|l is a particular case of Eq. (|10|l with G 
given by Eq. (|15|l at = for all k except for ai = 1/2. 
The soliton can be constructed by setting an arbitrary 
value for /_ (or /„) in the range [qm,qs] and finding /„ 
(or /_) from the quartic Eq. H36(l . Quantities qm and 
qs are the amplitudes of solitons centered between two 
lattice sites and on a lattice site, respectively. We have 
qm = V2uj, which does not depend on h and a. For 
/_ > qs Eq. does not have real solutions, i.e., qs 

corresponds to the magnitude of /_ for which the two 
distinct real roots of Eq. (|36(l in /„ merge into a multi- 
ple root. The arbitrariness in the choice of initial value of 
/_ (or /„) implies the absence of the Peierls-Nabarro po- 
tential and the possibility to place the soliton anywhere 
with respect to the lattice. 

The dynamical evolution of IV'nP for exact solitary 
wave solutions (constructed numerically as described 
above) in the classical DNLS of model / is shown in Fig. 
^ where we compare the solitons centered at (a) xq = 
and (b) xq = 1/2. The results are obtained by numerical 
integration of Eq. (|19|l for lattice spacing h = 0.4 and 
soliton frequency uj = 1. The stationary solution in (a) 
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FIG. 2: The dynamical evolution of IV'nP for exact solitary 
wave solutions placed non-symmetrically with respect to the 
lattice for (a) model II and (b) model /// with a — 0.2. 
In both cases the solitons do not radiate and they do not 
move to a higher-symmetry position since the Peierls-Nabarro 
potential is absent. Results for lattice spacing h — 0.4 and 
soliton frequency uj = 1. 

is stable while the inter-site centered one in (b) is un- 
stable; due to the presence of perturbations in the form 
of round-off errors, it spontaneously starts to alternate 
between the two nearest inter-site configurations passing 
through the stable on-site configuration. The pulse of 
Fig. nja) does not radiate while that of Fig. ^b) does. 

Similar computations have been carried out in models 
// and ///, and are reported in Fig. 13 also using h = 0.4 
and UJ ~ I. The solitary waves were found to be stable 
for any xq both in model // and model /// for negative 
and positive a in a vicinity of a = 0. For example, in 
Fig. 12 we show the results for (a) model // and (b) 
model /// with a = 0.2 and initial profiles placed non- 
symmetrically with respect to the lattice. In both cases 
the solitons do not radiate and they do not move to a 



higher-symmetry position. 



VII. CONCLUSIONS 

We have described a general and systematic method of 
constructing spatial discretizations of NLS-type models, 
whose stationary soliton solutions can be obtained from a 
two-point difference problem. In this setting, finding sta- 
tionary solutions becomes tantamount to solving simple 
nonlinear algebraic equations. We have also illustrated 
the connections of the resulting models with the inte- 
grable discretization of the NLS equation, of which they 
are a natural generalization for cubic nonlinearities (our 
construction was given for arbitrary polynomial nonlin- 
earities of a particular parity); furthermore, the differ- 
ences of such models from the standard discretization of 
the NLS equation often encountered in physical applica- 
tions have been highlighted, both in terms of the relevant 
dynamical (solitonic) behavior as well as in terms of the 
underlying conservation laws present in the various mod- 
els. 

It would be particularly interesting to further exam- 
ine such discretizations and their features, such as the 
stability of their solutions , and their travelling wave 
properties j3 IHl and potential integrability of spe- 
cial members within these families. Such studies are cur- 
rently in progress and will be reported in future publica- 
tions. 
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